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The Morpheus Lander is a vertical takeoff and landing test bed vehicle developed to 
demonstrate the system performance of the Guidance, Navigation and Control (GN&C) 
system capability for the integrated autonomous landing and hazard avoidance system 
hardware and software. The Morpheus flight control system design must be robust to 
various mission profiles. This paper presents a design methodology for employing numerical 
optimization to develop the Morpheus flight control system. The design objectives include 
attitude tracking accuracy and robust stability with respect to rigid body dynamics and 
propellant slosh. Under the assumption that the Morpheus time-varying dynamics and 
control system can be frozen over a short period of time, the flight controllers are designed 
to stabilize all selected frozen-time control systems in the presence of parametric 
uncertainty. Both control gains in the inner attitude control loop and guidance gains in the 
outer position control loop are designed to maximize the vehicle performance while ensuring 
robustness. The flight control system designs provided herein have been demonstrated to 
provide stable control systems in both Draper Ares Stability Analysis Tool (AS AT) and the 
NASA/JSC Trick-based Morpheus time domain simulation. 


Nomenclature 

(j) c = platform command angle (rad) 

(j) P = platform angle (rad) 

<f>c.g. = vehicle pitch angle about the center of gravity (rad) 
a = angle of attack (rad) 

Pe = engine gimbal angle (rad) 

Pc = engine gimbal command angle (rad) 

G gd = engine dynamics transfer function from gimbal command to gimbal angle (unitless) 
Z c g = displacement of vehicle c.g. normal to reference (ft) 

Z sj = sloshing fluid displacement of j th tank (ft) 

m sj = slosh mass in j th tank (slug) 

l sj = c.g. to the j th slosh-mass distance (ft) 

M = total vehicle mass (slug) 

Ci = c za qA(X cg -X cv )/I zz , aerodynamic acceleration coefficient (1/sec 2 ) 
c 2 = FX c g JI Z z , vehicle angular acceleration (1/sec 2 ) 

Izz = yaw vehicle moment of inertia (slug-ft 2 ) 

X c . g . = center of gravity measured from the gimbal (ft) 
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X c . p . = center of pressure measured from the gimbal (ft) 

S E = engine 1 st moment about gimbal (slug-ft) 

I E = engine 2 nd moment about gimbal (slug-ft 2 ) 

k 3 = ( F-D)/M , vehicle acceleration (ft/sec 2 ) 

k 4 = F/M, acceleration due to gimballing (ft/sec 2 ) 

ky = C Zc $A/M, acceleration due to aerodynamic forces (ft/sec 2 ) 

C Za = rigid-body aerodynamic side force coefficient slope (unitless) 
q = dynamic pressure (lbf/ft 2 ) 

A = reference area (ft 2 ) 

F = gimbaled engine thrust (lbf) 

V w = cross-wind velocity (ft/sec) 

V = vehicle velocity (ft/sec) 

Zpt = damping factor of the i th bending mode (unitless) 
ctyi = natural frequency of the i th bending mode (rad/sec) 

Mi = generalized mass of the i th bending mode (unitless) 

77, = bending displacement of the i h bending mode (slug 1/2 ft) 

Y pi = bending mode shape at gimbal (slug 172 ) 

V = bending mode slope at gimbal (rad-slug' 1/2 /ft) 

Y ft ’ = bending mode slope at platform location (rad-slug' 1/2 /ft) 

V = bending mode slope at rate gyro location (rad-slug' 1/2 /ft) 

(p 

Y sij = deflection of the i th bending mode at the j th slosh mass location (slug 172 ) 
C Zri i = wind bending force coefficient for the i th bending mode (lbf- slug' 172 /rad) 
n b = number of bending modes 

n s = number of slosh modes (determined by number of propellant tanks) 

= vehicle roll angle about the vehicle x-axis (rad) 

(j) r = vehicle roll rate about the vehicle x-axis (rad/s) 


I. Introduction 

The Morpheus Lander is a vertical takeoff and landing test bed vehicle 1 . The prototype vehicle was 
manufactured and assembled at NASA's Johnson Space Center (JSC). The Morpheus Lander is used to demonstrate 
the system performance of the Guidance, Navigation and Control (GN&C) system capability for the integrated 
autonomous landing and hazard avoidance system hardware and software. The Morpheus Lander consists of a 
throttleable main engine using liquid oxygen and liquid methane as propellant, four spherical aluminum tanks (2 for 
liquid methane and 2 for liquid oxygen) and the necessary avionics to perform powered descent and landing to 
provide realistic flight trajectories to demonstrate hazard detection and avoidance technology. The Morpheus flight 
control system utilizes Thrust Vector Control (TVC) of the main engine to control the vehicle Body pitch and yaw 
axes and fixed reaction control system thrusters to control the vehicle Body roll axis. The onboard GN&C sensor 
suite includes a Global Positioning System (GPS) receiver, International Space Station (ISS) version of Honeywell’s 
Space Integrated GPS/INS (SIGI), Litton LN-200 Inertial Measurement Unit (IMU) and an Acuity laser altimeter 1 . 
The Morpheus’ primary payload is the Autonomous Landing Hazard Avoidance Technology (ALHAT), which will 
scan and detect hazards at the targeted landing site and interface with Morpheus’ GN&C system to autonomously 
avoid landing on any hazards (e.g. boulder or crater). 

The Morpheus Lander has several characteristics, which must be taken into account in design. First, it is 
aerodynamically unstable because the vehicle’s center of pressure is above its center of gravity. Second, the 
vehicle’s mass, moment of inertia, aerodynamics and slosh characteristics vary greatly throughout the flight timeline. 
Third, the Morpheus flight control system must be robust to various mission profiles as shown in Figure 1 . The outer 
guidance loop interacts with the inner attitude control loop and may result in instability if excessive guidance gains 
are used. 

To achieve the desired performance and guaranteed robust stability for the Morpheus control system, the flight 
controller design can be divided into four phases: 
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(1) Optimal rigid controller designs, which include phase plane controllers for roll channel and Proportional- 
Integral -Derivative (PID) controllers for pitch and yaw channels. 

(2) Robust guidance gain design to guarantee robust stability and adequate performance. 

(3) Flex filter designs to guarantee robust stability and adequate performance. 

(4) High fidelity nonlinear simulation to verify stability and performance for entire Morpheus control system. 

These design phases are not necessarily performed in sequence. For example, guidance and controller designs can be 
combined into an iterative design process. 



A series of optimal guidance commands for each mission are designed for takeoff and ascent to altitude, hover 
and/or flight descent trajectory, and landing the Morpheus vehicle. The flight control system is designed to ensure 
stable response to guidance commands while minimizing trajectory deviations. The TVC system gimbals the main 
engine cone to minimize both pitch and yaw attitude errors. Slosh frequencies change over time during the entire 
flight. Flex filters in the flight control system are designed to minimize the slosh components in the error signals 
before they are sent to the attitude controller. 


Digital filter designs and application have been extensively investigated in the literature 2-4 . The most popular 
analog filter design approaches include Butterworth, Chebyshev and elliptic filters. Direct digital filter design 
methodologies are also available in the literature 3 . These methodologies have been used for open-loop system 
designs by shaping filters to meet the open-loop performance specifications in the frequency domain. Unfortunately, 
most of the closed-loop system stability/performance requirements cannot be directly mapped into open-loop system 
specifications. Recently, a robust controller design methodology was proposed using a numerical constrained 
optimization approach to maximize stability margins while meeting performance requirements 5 " 9 . This control 
design methodology has been successfully used to design a single robust control moment gyro (CMG) flex filter set 
for multiple International Space Station (ISS) stages 5 and robust flex filters and a PID controller for the Orbiter 
Repair Maneuver operation 7 . The same design methodology for employing numerical optimization was developed 
for the Ares I bending filters, PID gains and anti-drift/load relief controller. The design objectives include attitude 
tracking accuracy, load relief, anti-drift and robust stability with respect to rigid body dynamics, propellant slosh, 
and flex 10 . 


Under the assumption that the Morpheus time-varying dynamics and control system can be frozen over a short 
period of time, the bending filters are designed to stabilize all the selected frozen-time launch control systems in the 
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presence of parametric uncertainty. Due to the rigid-slosh interaction, the slosh filter design achieves stability by 
phase stabilizing the slosh modes. Similarly, the Roll Control System was designed to deal with disturbance torques 
due to aerodynamics and vehicle thrust vector and mass properties dispersions while providing steering torque for 
roll control commanded by guidance. These effects had to be considered in the design of the phase planes to provide 
the required attitude performance while conserving propellant. A design methodology for employing numerical 
optimization 5 ’ 10 was used for the Morpheus flight control system design. 

The Morpheus flight control system designs presented here have been demonstrated to provide stable control 
systems in Draper Ascent vehicle Stability Analysis Tool (AS AT) 11 and verified in the JSC Trick 12 Morpheus time 
domain simulation and flight tests. 


II. Morpheus Flight Control System 

The cross coupling effect between axes in the Morpheus Lander vehicle is small. This enables design and 
analysis to be performed in each axis separately. The Morpheus attitude control system consists of both gimbal 
control and fixed jet reaction control. The former is used for pitch and yaw control while the latter is utilized for roll 
control. Due to the symmetry in the pitch and yaw axes as shown in Figure 2, the TVC control system can be 
designed for either the pitch or yaw axis 13 . 



Figure 2. Morpheus Vehicle Coordinate Frame 


Lateral Control System 

The 2 DOF equations presented in their entirety in Reference 13-15 can be divided into rigid body dynamics, 
flex (bending) dynamics, and slosh dynamics. Figure 3 shows the sign conventions for the equations of motion with 
variables as defined in the List of Variables 13 . Note that the 2 DOF in this depiction include translation in the Y-axis 
and rotation in the Yaw direction. Rigid body dynamics includes rigid rotation about the center of gravity and the 
rigid body translation of the center of gravity. Specifically, 

&c.g. S —~C l OC—C 2 (^c.g^E ^e)Pe S )Zj/ (1) 

zz zz 7=1 

and 

z c .^ 2 = k 7 a + k 4 /3 E + k 3 0 cg + p E s 2 - A Yjn sj Z s / ( 2 ) 
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Figure 3. Spacecraft Sign Conventions 


Rigid body motion depends on forces and torques created by aerodynamic pressure, gimbal angle (including 
the “tail-wags-dog” effect), thrust, drag, and sloshing fuel masses. Angle of attack for the undeformed vehicle is a 
function of wind velocity (V w ), rigid attitude (</> c . g .) 9 and rigid velocity (Z cg s) 

a = (V w /V) + (j) cg - (Z cg s/V) (3) 

Bending motion depends on flexibility characteristics such as natural frequency, damping, and generalized mass in 
addition to the forces mentioned above for the rigid body dynamics. Specifically, 


(/ + 24 Al s + 4,)Mn = FY„p E + (44 - i e y'„)P e s 2 - 2>„ + c z , a 


Each bending mode is represented by a separate equation. The inclusion of more bending modes in a model leads to 
an increase in model fidelity. However, this also leads to an increase in the size and complexity of the system. The 
number of modes included in an analysis ultimately depends on the required model fidelity, the vehicle’s frequency 
spectrum, and the available computation tools. 


The dynamics of sloshing fuel is described by an equation for each fuel tank being considered in the analysis. 
Specifically, 

(s 2 + 24 + CO 2 )Z„ = - z, t s 2 + (/,/ + -f,Y„ A s 2 (5) 


Similar to bending modes, the inclusion of more slosh masses (up to the true number) yields a higher fidelity 
model. As shown in Equations 1, 2, and 4, the slosh tanks with greater slosh masses have a larger effect on overall 
vehicle motion. Equations 1 through 5 describe the overall yaw motion of the spacecraft and collectively form the 
plant dynamics. 


In addition to the plant dynamics, the control system includes a combination of other systems designed to 
regulate the spacecraft’s attitude. Sensors on the body of the vehicle determine the attitude and attitude rate 
according to the equations 

(6) 

1=1 

and 

n b t 

^g‘ v =4,,- v +Zv^ v (7) 

i=i 

respectively. It is apparent that sensor output only differs from plant output if bending dynamics are included in the 
model. Note that, when convenient, the sensor equations can be included as part of the plant’s output equations. 
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Sensor output is filtered in order to attenuate the higher frequency bending modes 10,13,14 . The filter output becomes 
the PID controller input. The PID controller outputs a gimbal command angle 

Pc ~ j* (^p -Qc ) ((ftp ~ (j) c ) + K D (j) RG (8) 

The controller gains ( K h K P , and K D ) are designed using gain scheduling, a process that involves optimizing a 
controller for specific design points throughout the time of flight and linearly interpolating between these points to 
obtain gain values for any flight time 10,13 . The command angle translates the gimbal angle via the dynamics of the 
actuator ( G gd ) 

Pe = G G dPc W 

The lateral attitude control system model, as shown in Figure 4 includes the Morpheus dynamics, guidance and 
controller modules. The dynamics modules consist of both rigid and slosh dynamics models 1 . The dynamics 
modules take in actual gimbal angles and output attitude and rate gyro measurements. The rate signals are filtered by 
the rate filters. The attitude commands are generated from the steering block, which takes in the weighted sum of 
position and velocity tracking errors. The attitude and rate error signals represent the difference between 
commanded and actual attitude and rate. The error signals are sent to the PID controller to generate the commanded 
gimbal angles, which drive the gimbal dynamics. 


Rotational Motion Control 



Figure 4. Morpheus Lateral Flight Control System Model 


Both position and velocity measurements are fed back in the outer loop. The weighted position and velocity 
errors are summed to generate estimated acceleration commands. The steering algorithm takes in the estimated 
acceleration commands to generate attitude commands. To improve attitude pointing accuracy, both lateral Y/Z 
accelerations are fed back in the outer loop. The acceleration filters are designed to update the trimmed attitude 
command from the sensed acceleration, as well as eliminate the high frequency sensor noise shown in Figure 4. 

Roll Control System 

The Morpheus flight control system uses phase plane autopilots for roll attitude control 1 . The roll phase plane 
autopilot produces a single command that takes on three states; 0 for drift, +1 for positive rotation and -1 for 
negative rotation. This command is sent to a jet selection logic that sets the appropriate fixed jet command(s) to “on” 
or “off.” 

Because the high-frequency flex dynamics possess the ability to make the system unstable, it is critical that a low 
pass filter be employed to attenuate high-frequency flex modes. The rigid and flex plant dynamics are integrated 


6 




with the flex filter and phase plane controller to create the roll dynamics and control system. The block diagram in 
Figure 5 outlines the dynamics and control components. 


Phase Plane 



Figure 5. Roll Control System Model 


The flex filter block consists of a set of low-pass filters to attenuate high-frequency noise while, at the same 
time, allowing low-frequency dynamics to feedback into the controller. A phase plane control system regulates 
attitude tracking and performance. Unlike the PID controller in the TVC control system, the phase plane controller is 
an inherently nonlinear system that necessitates the use of nonlinear techniques in order to predict the system’s 
behavior. The phase plane controller is “an idealized method of treating performance optimization for classes of 
minimum time and/or minimum fuel problems” 16 . The phase plane controller offers a unique method for attitude 
control while responding to the vehicle dynamics in the plane defined by state errors and state rate errors. The 
trajectories in the phase plane can be described through Equations 10 and 1 1 17 . 

(/> rl = </> r \ + CO At (10) 

fir! = //i + CX r At (11) 

In the expressions above co r is the roll angular velocity, a r is the roll angular acceleration and At is the thruster firing 
time for the phase plane controller. 

Figure 6 shows how trajectories in the phase plane operate. For example, consider the starting point in Figure 6. 
The system applies a continuous torque until it enters the drift channel. At this point, the thrusters discontinue their 
firing and the system’s attitude continues to increase because the system is in the upper half of the phase plane, 
which means it has a positive rate. The rate is constant as long as there is no firing because the system cannot 
accelerate. Once the system crosses the negative switch line into the negative firing region, the system undergoes a 
negative acceleration which places the system back into the non-firing region; however, the system’s attitude will 
continue to move towards the negative firing region until the system has been driven into the lower half of the phase 
plane. Once there, the system’s attitude will decrease in the non-firing region until it crosses the positive switch line 
at which point a positive firing will occur 17 . 
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The system will continue to oscillate around the origin of the phase plane in what is called a limit cycle. A 
common definition of a limit cycle is an oscillation of “fixed amplitude and fixed period without external 
excitation” 18 . In order to evaluate a nonlinear system such as a phase plane controlled system, it is necessary to 
transform the phase plane controller into a form where control techniques can be applied. This will be accomplished 
in a two-step process. First, only the attitude hold region will be evaluated, and second, the phase plane controller 
will be transformed into an equivalent system consisting of a PD controller and a nonlinear element consisting of a 
dead zone and an ideal relay. This procedure is illustrated in Figure 7. The two switching curves in Figure 6, which 
define the dead zone between the positive and negative firing regions, can be defined by the inequality: 

_ ^ ® RL (1^) 

In order to proceed from step one to step two in Figure 7, it is necessary to implement a PD controller. It is first 
necessary to recast Equation 12 in the form 

-1 <(/) r l 5 + ^) r ! CO^ <1 (13) 

Multiplying Equation 13 by the dead zone, 5, leads to: 

-8<4> r +(8l(o w yfr r <8 (14) 

This inequality provides the following values for K P and K D , which are the proportional and derivative gains 
respectively 17 . 

K p = 1 (15) 

K d =S/co rl (16) 

The K P and K D values calculated above are critical in creating a practical phase plane controller that can be used 
with control analysis techniques. Substituting these values into Equation 14 yields the result for the dead zone 17 . 

-8 < K p <j> r + K D (j) r < 8 (17) 

It has been demonstrated in Reference 19 that with the above approximation, both Pulse- Width Modulation (PWM) 
and absolute stability theory can be used to predict the stability margin of a phase plane control system. The 
applications of both techniques have been extended to analyze and design the Morpheus roll control system 19 . 

Original Phase Plane Control System 
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Figure 7. PD-Equivalent Phase Plane Development 


Morpheus Vertical Control 

The Morpheus vertical flight control system can be modeled as shown in Figure 8. The guidance block takes in 
the measured and commanded vertical position and velocity data, and then generates an acceleration command, 
which is the weighted sum of the position and velocity errors. Both commanded and sensed vertical acceleration 
data then send to the throttle control to generate a delta throttle command to compensate for the unmodeled 
acceleration force. A first order transfer function is used to represent actuator dynamics. 
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Cmd W*at 

Figure 8. Morpheus Vertical Flight Control System Model 


III. Constrained Optimization Approach for Morpheus Flight Control System Design 

It has been previously demonstrated in multiple space applications 511 that bending filters can be designed 
numerically using a constrained optimization framework. The design parameters are the coefficients of the bending 
filters. For example, if an n th order transfer function architecture is selected for both the attitude and rate filter, the 
total number of design parameters is 4n. 


F ( A = TT S + ^ X 4*i+2 X 4*i+l S + X 4*i+l 

" attitude A / f f 2 

1=0 S ~\~ ZXa, 


F ra,e ( S ) 


n 


S + ZX A *;, 0 Xa*;, 1 S + Xa, 


i = 0 S + 2x 4 : 


( 18 ) 


In addition to filter coefficients, PID and Guidance controller gains/filter coefficients are also included in the design 
parameters set. A set of feasible parameters must satisfy the following constraints: 

(Cl). The filter itself must be stable and minimal phase to guarantee stability and performance. 

(C2). The bandwidth of the bending filter should be greater than that of the PID controller to avoid rigid 
performance degradation. 

These constraints can be used to set the upper and lower bounds for the bending filter design. 

The primary objective of Ares -I control system design is to provide sufficient stability margins in the presence of 
various parameter uncertainties while maintaining adequate system response. Experience has shown that the 
following design objectives similar to Reference 13 are adequate: 

(01) . Both nominal and perturbed closed-loop Morpheus control system must be stable. 

(02) . At least 3 dB / 20 degree robust rigid gain/phase margin is required for perturbed control system. 

To ensure adequate response to guidance command, attitude step response specifications are imposed in the 
optimization problem. The additional constraints minimize the performance degradation due to the addition of the 
bending filters. These attitude pointing specifications, which include requirements for rise time, percent overshoot, 
settling time, and steady state errors, vary during the entire flight time. 
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The Morpheus control system must also ultimately demonstrate robustness to uncertainties in the plant. The goal 
is to design bending filters that are robust to uncertainty in slosh mass/damping/frequency, vehicle mass properties, 
engine mass properties, aerodynamics characteristics and transportation latency. Once design objectives and 
constraints are identified, the bending filter design task is ready to be cast as the following constrained optimization 
problem 


min fix) 

it. < 19 > 

g(x) ^ 0 

x,<x<x u 

The filter design criteria (Cl) and (C2) can be formulated as inequality constraints. These inequality constraints 
can be also cast as objectives f(x) in the above multi-objective constrained optimization problem. In general, these 
objectives are competing with each other. For example, maximizing gain margins usually diminishes phase margins. 
Therefore, there is no unique solution to this problem. To address this, Pareto optimality 20 must be applied to 
characterize the objectives. This is accomplished with a weighted sum strategy, which converts the multi-objective 
problem into a single objective optimization problem. 

IV. Analysis of Morpheus Flight Control System Design 

The Morpheus flight control system consists of Thrust Vector Controllers (TVCs) using identical pitch/yaw 
autopilot architectures and a phase plane autopilot for roll control. 

The thrust vector control law is PID with gains based on mass table lookup to accommodate changing flight 
conditions. The PID gains were designed for rigid body performance, but also to provide as much frequency 
separation as possible between the rigid body and slosh mode frequencies. The flex filters were designed to phase 
stabilize the slosh modes. Slosh dynamics are also damped passively using mechanical baffles in the tanks, and 
filtering in the control system is used to supply adequate slosh gain/phase margin for phase stabilized slosh modes 
when passive damping is inadequate for slosh gain stabilization. 

The robustness of the Morpheus flight control system designs are verified in the linear frequency analysis using 
Draper Ascent vehicle Stability Analysis Tool (AS AT) for all dispersed designated trajectories. An AS AT utility 
imports each nominal datum and identifies the correct slosh mass property for each of a series of frozen time points. 
The analysis time vector is the same for each trajectory and is chosen to provide snapshots of flight with sufficient 
control authority (thrust magnitude) such that the stability analysis is valid. The margin calculation tool determines 
the characteristics of the response and tracks both the worst-case margin values and the seed identifiers that created 
the worst-case systems such that they can be later replicated. The following margins and characteristics are 
identified if they exist: 

Rigid gain margin: minimum gain margin at rigid phase crossover frequencies. 

Rigid phase margin: minimum phase margin at rigid gain crossover frequency. 

Rigid slosh margin: minimum gain margin at slosh phase crossover frequencies. 

The relevant characteristics are illustrated on the Nichols response shown in Figure 10. The robust frequency margin 
requirements for the current design are listed in (02). 

An example of a stability margin evaluation of the nominal and dispersed control system using ASAT analysis 
are summarized for both the vertical control system (Figure 9(a)) and the lateral control system (Figure 9(b)). The 
design goal is to ensure that the Morpheus flight control system possess sufficiently robust stability gain and phase 
margins for the entire flight region of multiple planned missions. 

After the flight control system is designed and analyzed in the frequency domain, a high fidelity nonlinear Trick 
simulation is used to verify stability and performance for the Morpheus control system. A Monte Carlo analysis is 
performed by perturbing system parameters such as mass property, center of mass location, thruster misalignment, 
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as well as biased wind to evaluate the performance and robustness of the Morpheus flight control system. A sample 
Monte Carlo Analysis result is illustrated in Figure 1 1 . 
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Figure 9. Morpheus ASAT Stability Analysis Example 
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Monte Carlo Truth Trajectories 
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Figure 11. Morpheus Trick Monte Carlo Analysis Example 


V. Summary 

The Morpheus Lander is a vertical takeoff and landing test bed vehicle developed to demonstrate the system 
performance of the GN&C system capability for the integrated autonomous landing and hazard avoidance system 
hardware and software. This paper presents a design methodology employing numerical optimization to develop the 
Morpheus flight control system. The flight control system design methodology was based on a numerical 
constrained optimization approach to maximize performance while meeting stability margins requirements. The 
flight control system designs provided herein have been demonstrated to provide stable flight control system in 
Draper AS AT, NASA/JSC Trick Morpheus time domain simulation, and flight tests. 
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